Cavity assisted quasiparticle damping in a Bose-Einstein condensate 
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We consider an atomic Bose-Einstein condensate held within an optical cavity and interacting 
with laser fields. We show how the interaction of the cavity mode with the condensate can cause 
energy due to excitations to be coupled to a lossy cavity mode, which then decays, thus damping 
the condensate, how to choose parameters for damping specific excitations, and how to target a 
range of different excitations to potentially produce extremely cold condensates. 
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The breakthrough success, largely due to improved 
cooling techniques H, of producing atomic Bose-Einstein 
condensates (BEC) [pi, has created an exciting new area 
of investigative study [|[ [Ij. Within the field of cavity 
quantum electrodynamics, one is able to precisely con- 
trol the interaction of the electromagnetic field with an 
atom Q; the trapping and cooling of a single atom in 
an optical cavity has been successfully experimentally 
demonstrated 0. Here we combine |7j the use a cavity 
to cool atoms S with the ever present goal of producing 
colder condensates || , particularly relevant to control- 
ling quasiparticle excitations [|0, 11, l3] produced in the 
course of other experimental investigations [fl3|| . 

We consider a Fabry-Perot optical cavity with decay 
rate k, and a pump laser, with a BEC held in a har- 
monic potential of frequency u) inside the cavity, and 
driven by a separate laser O, as shown in Fig. [l]. We 
generally consider only one spatial dimension, consid- 
ering the radial motion to be frozen out by a tight 
harmonic potential of frequency to r (cigar-shaped con- 
figuration Jl4| ). The laser-cavity superposition field is 
E(x,t) — fl(x)e~ luJLt + gacos(kx) where a is the cav- 
ity mode annihilation operator; interacting with this are 
the BEC atoms, considered to have two levels \g) and 
|e), with transition frequency ujq. We initially consider 
a single particle, and ignore the degrees of freedom due 
to the free cavity field and the atomic motion. In a ro- 
tating frame (U = exp[— i(a ee + a^a)ujLt}), the Hamil- 
tonian describing the internal atomic level dynamics is 




FIG. 1: Proposed configuration: a BEC held in a cavity, 
driven by a laser f2 and a pump laser, and losing cavity pho- 
tons with a decay rate k. 
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,+E(x)a eg +E^(x)a ge ], where cr eg = \e)(g\, 
w cx , and E(x) = £l(x) + ga cos(fcx). Assum- 
E, A a 3> 7 where 7 is the atomic sponta- 
neous emission rate, and |A|, where A = lj c — uj cx , by far 
exceeding the frequency scales governing atomic motion 
and the cavity dynamics, we may adiabatically eliminate 
\e) to derive a closed time-evolution equation for the \g) 
wavefunction component [ jTsL using this to reconstruct 
an effective Hamiltonian: 
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Assuming g <^.Q, and (a' a) < 1, we can neglect the terms 
oc g 2 , and as 7 <C A a , we ultimately consider only the 
Hermitian terms. The light field thus provides: a term 
oc |r2(z)| 2 , which, assuming fl(x) oc e lkhX , contributes 
only a shift; and a contribution linear in a and a\ which 
acts like a driving field for the cavity mode. It thereby 
implements a mechanism to transfer energy to the cavity 
which in turn will shed this energy to the environment 
by means of cavity damping. The effective single particle 
Hamiltonian (including the free dynamics of the particle 
motion and the cavity) is now 
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where f(x) = — Q(x)g cos(fcx) A a /(A 2 + 7 2 ), and m is 
the particle mass. We now consider a BEC-cavity sys- 
tem. The one-dimensional atomic interaction potential 
= u5(x — y), where u — 2huj r a s and a s is the s-wave 
scattering length H. Thus, in second quantized form: 



H cB = f dx¥{x)\n[f{x)a + /*(s)o+] 
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where $(x) and W(x) are atomic field operators, and 
the a p terms are provided by the pump laser. We treat 
^>(x) and a semiclassically, replacing them with the scalar 
quantities ip(x) and a, respectively. It is convenient to 
rescale to dimensionless harmonic units {% = m = u) = 
1); this scaling is assumed from now on. The result- 
ing equations of motion are: a Gross-Pitaevskii-like (GP) 
equation 

i<p{x) = {H GP + [f(x)a + f*(x)a*]}tp(x), (4) 

where Hqp = —d 2 /2dx 2 + x 2 /2 + v\tp(x)\ 2 is the un- 
perturbed GP "Hamiltonian," with v — uNy/m/hcu/h, 
f(x) = \/~Nf(x)/oj, and a — a/VA, where N is the par- 
ticle number; and the semiclassical cavity mode equation 
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We set a p to be = J dx\tpo(x)\ 2 f* (x), so that when ip(x) 
approaches ipo(x) (the ground state of Hqp), a will sim- 
ply decay, without feeding into Eq. (Q). The cooling 
mechanism thus switches off upon reaching the desired 
steady state <p{x) = tfo{x) and a = 0. 

We now consider linearized perturbations 8(p(x), 8a, 
around the steady state. Linearizing Eq. (kj) produces 
Bogoliubov-like equations E3j. Note in Eq. S) for ipa(x) 
no distinction is made between different global phases; 
in calculations the final phase of fo(x) is determined by 
the chosen initial conditions, and is not relevant. We 
consider only perturbations orthogonal to ipo (x) , and can 
implicitly assume a U(l) gauge transformation JlTJ] such 
that tpo(x) is always real. Having defined 8tp(x) as or- 
thogonal to <fio(x), i.e. S(p(x) — j dyQ(x,y)S(p(y) where 
Q(x,y) = 8(x — y) — ipo(x)ip (y), we can state that 

-/*«*»>( It)) 

+ J dyQ(x,y)[f(y)6a + r(y)Sa*} ( ) , 

where £(x,y) = J f dzdwQjx, z)H Bos (z, w)Q(w, y) 
in terms of the usual Bogoliubov Hamiltonian JlTJ] : 

_ / H GP + v\tp (x)\ 2 
~ ' ~v V *(x) 2 
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where fx is the ground state chemical potential; and 
It is convenient to expand (8tp(x), 8ip*(x)) as: 
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where {(uk(x),Vk{x)), (vk(x),Uk{x))} are eigenstates of 
L{x,y), with eigenfrequencies ±Ufc (these may be deter- 
mined numerically by diagonalizing C(x,y) @). As (p Q 
is assumed real, these are also real. We similarly expand: 



dyQ{x,y)f(y) 



k 



<po(y) 
-<po(y) 

u k (x) 
v k (x) 



v k (x) 
u k (x) 



where the coefficients can obviously be defined by 



Xk = dx[u k (x) +v k (x)]f(x)(p (x) 
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Note that treating the evolution of perturbations or- 
thogonal to (po(x\ is equivalent to a number-conserving 
formalism pl[ The Uk{x) and Vk(x) are orthog- 

onal to ipo(x), and the {(u k (x), v k (x)), (v k (x), u k (x))} 
are used as a convenient time- independent basis. All 
time dependence in Eq. (||) is thus in the Cfe, Q coef- 
ficients, in contrast to Refs. [jll], [l^, We now trans- 
form Eq. (||) to an interaction picture. We thus set 
(6<p(x),6<p*(x)) = e lCt (8<p(x),8tp*(x)) and 8a = e zAt (5a. 
Substituting Eqs. (§) and @ into an appropriately 
transformed Eq. (^|), and then making the integration 
J dxui(x)8(p(x) — vi{x)8(p (x), we end up with 

if, = texie 4(w! - A)t + te*xre l(wi+A)t , (12) 



where Q = e luJlt Xi- From Eqs. (||) and (10), the linearized 
equation of motion for 8a, after adiabatic elimination 
(assuming k \PkXk\)j produces 
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We assume that due to our choice of f(x), \xi \ dominates 
|Xfc^;|j neglecting terms where k ^ I, and choose A = 
to match it. Applying the rotating wave approximation 
(RWA) to Eqs. @ and @, we neglect terms involv- 
ing f* and <5a* (2cj, should therefore describe the fastest 
timescale). Combining the resulting expressions, we end 
up with a simple damping equation: 
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As fl(x) tx e^ 1 - 31 , we set f(x) = gQ[e*""> x " + e" 
where fc cx = k ± fc^, and A: is the wavenumber of the 
cavity mode. Taking only the e lk °* x term, we see in 
Fig. ||b) that each \xi\ 2 peaks for some k cx at a point 
where \xi\ 2 > \xk^i\ 2 - Assuming cjj is chosen to match, 
for sufficiently large k the e l ( k <=^ 2k ) x term can be ignored, 
as those modes to which it couples most strongly (i.e. 
large xi) can i by the RWA, be neglected. From now on we 
consider f(x) = g e lkcxX . As e lkcxX <f (x) is equivalent to 
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FIG. 2: a) Plots of |xi| 2 f° r ' =1-15 as functions of k ex , where 
Xi is as defined in Eq. (ll]), and fix) = e lkcxX . The vertical 
lines mark the numerica ly determined optimal values for fc e x- 
b) Plots of: fcox = V^ioi, i.e. approximately optimal value of 
fcox (plusses), and numerically determined optimal vales of 
fccx (circles), against u)i for I =1-25. 
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a momentum-kicked ground state with mean additional 
energy = fc 2 x /2, the inner product of f(x) (ipo (x) , —<po(x)) 
with eigenstates of C(x,y) should be large with those 
eigenstates of similar energy; one would thus simplisti- 
cally expect |xz| 2 to peak when fc 2 x /2 = uj\. In Fig. ||b) 
we see that this prediction is a little crude for small u>l, 
although it converges to the true values for higher fre- 
quencies, as expected. Also correctly predicted is that 
the optimal fc ex values tend to converge for large /. From 
Fig. ||a) we see that the |xz±i| 2 (at least) can be signif- 
icant compared to \xi\ 2 at the optimal value of fc cx , but 
the fact that these terms are both smaller and rotating 
means that the /th term will still dominate. For large I, 
as the peak values of k cx converge and the peaks become 
less well resolved, Eq. ([hi]) will become less reliable. 

Numerically we monitor the damping via the energy 
E = fdxtp*(x)[-d 2 /dx 2 + x 2 + v\<p(x)\ 2 }<p(x)/2. We 
define Eq equivalently, with <p(x) = tpo(x). Within the 
linearized approximation E = Eo+J2k w fc|Cfc| 2 > an d thus, 
if initially only Q ^ 0, from Eq. (Q) we have 



(15) 



where SE = E — Eq. The energy damping timescale 
is thus = 2g 2 A/n, where A = \xi\ 2 for g = 1. In 
Fig. ^ we compare Eq. ( |l5| ) with numerical integra- 
tions of Eqs. (||) and (Eh. The initial conditions are 
f(x) oc (fo(x) + 0A[ui(x) +vi(x)], a = 0; v is always = 10. 
In Fig. ||a) we see qualitative agreement of the analyti- 
cal estimate given by Eq. (15) with the numerical results 
for different go, which, in view of the number of approx- 
imations made, is remarkably good. In Fig. [|b) we see 
the convergence of the position density p(x) — \ip(x)\ 2 
to |(^o(a;)| 2 - In Fig. ||c) we see good qualitative agree- 
ment in the damping rates for g$ and k = 1, 10 with 
the analytical prediction (which is the same for each). In 
Fig. ||d), where I — 1 rather than 12, but the parame- 
ters are otherwise the same, we see heating for the k = 10 
case; u>\ = 1 is dominated by k, and the RWA is therefore 



FIG. 3: a) Plots of SE/Eo for different g% against time, 
where / — 12, n — 10. Solid lines, numerical calculations; 
dotted lines, prediction of Eq. jl^). b) Time evolution of 
p{x) = \tp(x)\ 2 (I = 12, k, go = 10). c) Plots of 5E/E 
against time, where / = 12. Solid line, k, g% = 10; dashed 
line, k, g% = 1; dotted line, prediction of Eq. ( |l5| ) for both 
cases, d) Plots of SE/Eo against time, where I — 1. Solid line, 
K, go = 10; dashed line, k, g% = 1; dotted line, prediction of 
Eq. ©. 



not valid. Most significant about equations ( |l4| ) and jl5| ) 
is that the procedure for deriving them provides a recipe 
for determining parameters optimal for the damping of 
specific excitations. 

In the case of a finite temperature BEC, where there 
is a range of populated quasiparticle (QP) excitations, 
these can be individually targeted, resulting in cooling of 
the BEC. There is a well known correspondence between 
linearized perturbations of <p{x) and QP excitations |ll], 
|i"2] |. Assuming ajdo ~ N, where ao is the annihilation 
operator for a particle in state (pa(x) and N is the total 
particle number operator, then 



*0) = ao 



(p (x) H = V" b k u k {x) + btv k (x) 



, (16) 



where the b k , b k are QP creation and annihilation op- 
erators, respectively. If the state of the system is as- 
sumed thermal, then (H^fc) = [exp(ojfc/r) — 1] _1 and 
(bk) = (by = 0j where r = fc^T/Titj. Semiclassically 
one can regard these expectation values as describing the 
statistics of Gaussian random variables (3k, with mean 
and variance (btbk) flql - One can make use of the 
linearized perturbation-QP correspondence by taking a 
random initial condition: 



ip(x) = (po(x) + -= Y]/3 fc Wfe(x) +/3fcU fc (x), 
v Jy 



(17) 
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FIG. 4: a) Damping of states perturbed by a single excita- 
tion, for ft, (?q = 1. Solid lines, numerical calculation; dotted 
lines, analytical estimates, b) Plots of SE/Eo against time: 
solid line, r = 10; dashed line, same initial condition, but 
with a cutoff for I > 25; dotted line, r = 5. c) White, (bjbi) 
for r = 10. Data corresponding to the solid line of b): grey, 



\Pi\ 2 for t = 0; black, |/3;j 2 for t = 100. 



with a = 0, and simulating the cooling process numer- 
ically with Eqs. ([|), (g). As we must target a range of 
excitations, we employ a procedure so that whenever t is 
a multiple of 4, k cx and A are changed (by adjusting the 
laser fi), initially targeting I = 25 (as described above), 
and then through each of I =24-1. We choose <?q, n = 1 
(to ensure cooling takes place for small u>i), N = 1000, 
and v = 10. In Fig. |Ja) we show damping of states of 
the form tp(x) oc <po(x) + 0.1[u/(x) + vi (x)], for a range of 
relevant I. For I = 24 there are deviations from Eq. (15), 
as expected; significant damping nevertheless takes place 
for each I. In Fig. ||b) we observe significant damping 
for initial states where r = 5, 10. Note that the energy 
separation between the r = 10 state, and the "cutoff" 
state [k = 1-25 in Eq. (|l7|), rather than 1-100], remains 
relatively static, implying that damping is taking place 
in a targeted manner on the low-lying excitations. Note 
also from Fig. ||c) that although the relative population 
for k > 25 is miniscule, the relative energy contribution is 
significant. Although the cutoff state was propagated for 
the sake of comparison, it could conceivably apply to a 
physical situation, such as occurs in evaporative cooling 
101. In Fig. |c) we see clear supression of the \Pi\ 2 , where 



(3i is determined by \A/V J dx<p* (x)ui(x) — ip(x)vi(x). By 
our semiclassical analogy, it follows that QP populations 
will similarly decrease, potentially providing a method to 
obtain extremely cold condensates. Note that some spon- 
taneous emission is inevitable, and the resulting heating 
must be slower than the overall cooling rate. Details on 
the effect of spontaneous emission on BEC temperature 



are in Ref. @. 

We have presented a cavity-laser-BEC configuration 
whereby, through astute choice of parameters, specific 
excitations of the BEC can be rapidly damped. We have 
shown the derivation of a simple equation which explains 
all major features of the damping process, and provides 
a recipe for determining optimal parameters. We have 
demonstrated how this damping procedure can in prinici- 
ple be used to produce extremely cold condensates. 
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